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In the last decade the CDF experiment at the Tevatron clearly demonstrated that it is possible 
to study extensively heavy flavour physics in hadron collisions and achieve remarkable results, 
competitive and complementary to B-factories. In this paper we report on the indirect searches 
for physics beyond the standard model via measurements of rare fo-hadron decays. The final 
limits, based on the analysis of the full CDF data set, on the branching fraction of the B°,s decay 
into a pair of muons are presented and discussed. Moreover we review the latest measurements, 
with 6.8 fb _1 of collected data, of the total and differential branching fractions and angular 
observables of rare fo-hadron decays proceeding via the flavour-changing neutral-current process 
b sjj. + Li~. 
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1. Introduction 

Rare flavour-changing neutral-current (FCNC) 2?9 * — > and b — > spL + pL~ decays are con- 

sidered among the most promising probes of the standard model (SM) and its extensions. The 
precise measurement of several observables (total and differential branching ratios, angular dis- 
tributions of the decay products) of these decays might provide interesting clues for new physics 
(NP) phenomena, if any sizable deviation from the SM predictions is observed. In this paper we 
review the current status of these indirect searches at the Collider Detector at Fermilab (CDF II), 
which reached a sensitivity very close to SM predictions after a decade of Tevatron leadership in 
the exploration of B® dynamics. The Tevatron pp collider, whose operations ended in October 
2011 after 20 years of operation, provided excellent opportunities to study B physics. CDF II is 
a multipurpose detector, consisting of a central charged particle tracking system, surrounded by 
calorimeters and muon chambers. It collected a final data set corresponding to about 10 fb _1 of 
integrated luminosity. 

In the search for rare B decays, the experimental challenge is to reject a huge background while 
keeping the signal efficiency high. A dedicated dimuon trigger has been used to select events with 
a pair of muons in the final state in the pseudorapidity region |tj | < 1 . 1 . In the reconstruction and 
analysis of B hadron decays, CDF takes advantage of an excellent transverse momentum resolution, 
^7 = 0.07% pr (GeV/c), which implies a resolution on the invariant dimuon mass of 24 MeV in 
the B9s — > £i + £t~ decay, a vertex resolution of about 30 /im in the transverse plane, and particle 
identification (PID) capability, based on multiple measurements of the ionization per unit of path 
length (dE / dx) in the drift chambers. 

2. Search for fi9 n -»■ pL + pr 

B9 s — > £i + jU~ decays are mediated by FCNC and thus forbidden at first order in the SM. More- 
over they are further suppressed by helicity factors (m^/niB) 2 in the final dimuon state. They can 
only occur at second order through penguin and box diagrams. The SM predicts very low rates for 
these processes: m(B° s -> jU + aO = (3.2 ±0.2) x 10~ 9 and ^(5° -> /i+ju - ) = (1-0 ±0.1) x 10~ 10 
However, a wide variety of beyond the standard model (BSM) theories predict enhancement 
of their branching ratios by several order of magnitudes, making these decays one of the most 
sensitive probes in indirect searches for NP. 

In 201 1, CDF observed an intriguing ~2.5a excess over background in B° s — > jU + /i~ using 7 
fb _1 of data J^. Though it was compatible with other experimental results (LHCb (Q), CMS 
and the SM prediction, it could be interpreted as the first indication of a signal and allowed CDF 
to set a two-sided bound on the rate S3(jf s — > jU + /0 = (1.8+Q.j,) x 10~ 8 . To further investigate 
the nature of the excess, we repeated the analysis unchanged using the whole Run II data set, 
corresponding to 9.7 fb _1 of integrated luminosity, about 30% more data with respect to the 2011 
analysis. Here we report on the final results of this search 

The baseline selection requires high quality muon candidates with opposite charge, transverse 
momentum pj > 2 GeV/c, and a dimuon invariant mass in the range 4.669-5.969 GeV/c 2 . The 
muon pairs are constrained to originate from a common, well-measured reconstructed decay point. 
A likelihood-based muon identification method, is used to suppress contributions from hadrons 
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Figure 1: Dimuon mass distributions for the (a) B° — > and (b) B® — > signal region in the 

eight NN bins. The observed data (points) are compared to the total background expectation (light gray 
histogram). The hatched region is the total uncertainty on the background expectation. In (b) the dark gray 
histogram represents the SM signal expectation enhanced by a factor 4. 1 . 
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misidentified as muons. The branching ratios of B ^ — > pL + pL~ are measured by normalizing to 
a sample of 40225± 267 B + — > J/\j/(^- pL + pL~)K + candidates, selected with the same baseline 
requirements. A Neural Network (NN) classifier is used to improve signal over background sep- 
aration. Fourteen variables are used to construct the NN discriminant that ranges between and 
1. The six most discriminating variables include the 3D opening angle between the dimuon mo- 
mentum and the displacement vector between the primary and secondary vertex; the isolation / 1 
of the candidate ; the muon and B ^ impact parameters; the S9 decay length significance; the 
vertex-fit % 2 . The NN from the 2011 analysis was used with the same training. The final search 
region in the dimuon invariant mass has a half width of about 60 MeV corresponding to 2.5 times 
the dimuon mass resolution. The NN was validated with signal and background events. Careful 
checks for possible mass-biases of the NN output and overtraining show no anomalies. 

Extensive and detailed background estimates and checks have been performed. Background 
is due to both combinatorial and peaking contributions in the signal region. Combinatorial back- 
ground is estimated by fitting the sidebands to linear functions, after blinding the signal region in 
the dimuon mass distribution. The peaking background is due to B — > h + h'~ decays where the 
hadrons (h, h' stand for % or K) are misidentified as muons. This has been estimated from both MC 
and data. The misidentification probability is parametrized as a function of the track transverse mo- 
mentum using D*-tagged D° — > n + K~ events. It turned out that the peaking background is about 
10% of the combinatorial background in — > jU + jU~ and about 50% of the total background 
in the B° — > /i + jlt~ channel. Background estimates have been cross-checked using independent 
background-dominated control samples, in which the muons have the same measured charge or 
the reconstructed dimuon candidate lifetime is negative. No significant discrepancies between the 
expected and observed number of events in the control samples have been found. 
The data are divided into 8 bins of the NN discriminant to exploit the improved background sup- 
pression at high NN values, and five bins of mass in the search region. In the B search region data 
are consistent with the background prediction (Fig. 1(a)) and yield the limit of ^(5° — > pL + pL~) < 



3.8(4.6) x 10~ 9 at 90% (95%) C.L.. The significance of the background-only hypothesis expressed 
as ap-value, estimated from an ensemble of background-only pseudo-experiments, is 41%. 

In the B® search region, a moderate excess in the highest NN bins (>0.97) is observed (Fig. l(b] ) 
The p-value for the background-only hypothesis is 0.94%. We also produce an ensemble of simu- 
lated experiments that includes a 5° — > /i + jU~ contribution at the expected SM branching fraction 
which yields ap-value of 7.1%. With respect to the previous CDF result with 7 fb _1 of data [pj|, the 
excess in the third significant NN bin (0.97-0.987) softened, as expected for a statistical fluctuation. 
Though the 201 1 hint of signal is not reinforced by the new data, it is still present and remains >2a 
significant over background. Assuming the observed excess in the B® region is due to signal, CDF 
finds 3§(B° S -> jU + aO = (I.3+07) x 10~ 8 , which is still compatible with both the SM expectation 
and the latest constraints from LHC experiments [0, 

3. b -»■ decays 

Rare decays of bottom hadrons mediated by the FCNC process b — > s/j. + /j.^ are suppressed 

'/ = IPj-^l/CE,! Pt + IPr^D' w here p 1 ^ is the momentum of the dimuon pair; the sum is over all tracks with 
VWF-R^tF < I". At) and A<j> are the relative azimuthal angle and pseudorapidity of track i with respect to . 
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CDF Run II Preliminary L=6.8fb 
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Figure 2: Invariant mass of (a) B + -> K + fi+^-, (b) B° -> K*°pL + jX-, (c) 5° ->■ 0£i + jU~, (d) fi+ -> 
K* + ii + n~, (e) B° -» K s n + ii~, (f) A° -» An + n~ with fit results overlaid. 



at tree level in the SM and must occur through higher-order loop amplitudes. Their expected 
branching ratios are of the order of 10~ 6 . Because of their clean experimental signature and the 
reliable theoretical predictions for their rates, these are excellent channels for NP searches. 

CDF has studied the FCNC — > hji + }i~ decays (where Hi, and h indicate hadrons containing 
a b and s quark, respectively) listed in Table [j], using 6.8 fb _1 of data collected with the dimuon 
trigger Candidates for each decay have been selected by standard kinematics cuts and a NN 
optimized for best sensitivity. Signal yields are obtained by an unbinned maximum log-likelihood 
fit to the ^-hadron mass distributions (Fig. ||]), modelling the signal peak with a Gaussian, and the 
combinatorial background with a linear function. To cancel dominant systematic uncertainties, the 
branching ratio of each rare decay Ht, — > hjx + n~ is measured relative to the corresponding resonant 
channel H\, — > J/\j/h, used as a normalization and a cross-check of the whole analysis. The results 
of the total branching ratios are reported in Table [j] and include the first observation of the baryonic 
FCNC decay ->■ AjU + /i _ , and the first measurement of the B + -4 K* + ii + [i~ and B° -> K s 
decays at a hadron collider. 

Rich information about the b — > spL + pL~ dynamics can be obtained by precise measurements 
of the differential branching ratio as a function of q 2 = mlLc 2 and the angular distributions of 
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Decay mode ^(KT 6 ) 



B+ -^K+ix+^i- 
B° s ^<l>n+n- 

Ag^A/x+/i- 



0.46 ±0.04 ±0.02 
1.02 ±0.10 ±0.06 
1.47 ±0.24 ±0.46 
0.95 ±0.32 ±0.08 
0.32 ±0.10 ±0.02 
1.73 ±0.42 ±0.55 



Table 1: Branching ratio of the Hb —> decays measured by CDF. The first quoted uncertainty is 

statistical, the second is systematic. 
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Figure 3: Differential branching ratios of (a) B+ -> K + n+n~, (b) B° -)• K*°n + jj,-, (c) B° -)• , (d) 

B + -> K* + n + n~, (e) B° -> H + fi~, (f) AjJ -> AjU + jU~. The points are the fit result from data. The solid 
curves are the SM expectation [|ll] |T^, |l3|, [l4|]. The dashed line in (f) is the SM prediction normalized to our 
total branching ratio measurement. The hatched regions are the charmonium veto regions. 



the decay products. The differential branching ratios with respect to q 2 have been measured by 
dividing the signal region into six bins in q 2 and fitting the signal yield in each bin. In each fit, 
the mean of the Hy mass and the background slope were fixed to the value from the global fit, so 
that only the signal fraction was allowed to vary in the fit. The results are shown in Fig. ||. For 
B° s — > §n + ii~ and A° — >• A.jX + }X~ these are the first such measurements. At present no significant 
discrepancy from SM prediction is found. 

The angular distributions of the combined B° — > K* Q n + n~ and B + — > K* + pL + }l~ decays have 
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Figure 4: Measurements of angular observables (a) Afb, (b) Fl, (c) A^, , and (d) A,,„ as a function of dimuon 
mass squared q 2 in the combined decay mode B — > K* ji + n~. The points are the fit results from data. The 
solid and dotted curves represent expectations from the SM and a particular BSM scenario, respectively. 



been measured and parametrized to four angular observables: the muon forward-backward asym- 

(2) 

metry Aps, the K* longitudinal polarization fraction Fl, the transverse polarization asymmetry A T , 
the time-reversal-odd charge-and-parity asymmetry A\ m , defined in [|l^,[T^]. A j and A, m have been 
measured for the first time by CDF [|To|]. The results for these observables, shown in Fig. |], are 
among the most precise to date and consistent with SM predictions and other experiments, but still 
statistically limited in providing stringent tests on various BSM models. 



4. Conclusion 

We have summarized the recent updates on the searches of rare ft-hadron decays at CDF. The 
intriguing excess in fi° — > jU + /i _ reported in 2011 is confirmed with the full data set, though its 
significance is softened to the level of 2a over background. The measured B° — > jU + jU~ branching 
ratio is still compatible with the SM expectation and recent combined results from LHC experi- 
ments. 

The dynamics of several rare decays mediated by the FCNC process b — > sll^/j.^ has been 
studied in detail extending the reach to new angular observables. The A® — > A/i + jU~ decay has 
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been observed for the first time. Analyses of the b — > sjl fl decays are still in progress and may 
yield interesting results in the near future. 
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